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Abstract 

The recently measured direct CP asymmetries in the processes D^ — )• vr^vr" and 
fS| D^ — )• K^K~ show a significant deviation from the naive Standard Model expectation. 

^^ Using a general parameterization of the decay amplitudes, we show that the measured 

^vg branching ratios imply large SU{'i) breaking and large violations of the naive 1/A^c 

^ counting. Furthermore, rescattering constrains the / = amplitudes in the vrvr and 
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KK channels. Combining all this information, we show that, with present errors, the 
observed asymmetries are marginally compatible with the Standard Model. Improving 



^ the experimental accuracy could lead to an indirect signal of new physics. 

^ 1 Introduction 

^ The LHCb collaboration reported an interesting result based on 0.62 fb~^ of data at the 

^ Hadron Collider Physics Symposium 2011 [I]: 

(N AAcp = Acp{K+K-) - Acp(7r+7r-) = [-0.82 ± 0.21 (stat.) ± 0.11 (sys.)] %, (1) 

> 

which deviates from zero at 3.5cr level. Note that the effects from indirect CP violation cancel 
to a large extent in the sum, and a non-vanishing /S.Aqy' originates from the difference of the 
direct CP asymmetries, as explained below. In the above expression, the time-integrated CP 
asymmetry Acp(/) may be written as follows due to the slow mixing of neutral D mesons [2]: 



^cp(/) - ^g:;i;g:;-i - 4?(/) ^ <^ /"* ^ ^/(*) - ^-m/) + if « 



CP ' 



X (2) 

where Df{t) is the observed distribution of proper decay time and T£,o is the lifetime of the 
neutral D mesons. The indirect CP- violation parameter is given in terms of the parameters 
X = Aitid/Td, y = ATd/{2Td), \q/p\ and = arg(g/p) pi: 
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where r^cp = +1 is the CP parity of the final state considered here and (f) is the CP- violating 
phase. The HFAG average of the indirect CP asymmetry is A^ = (0.123 ± 0.248) % (i}]?]. In 
addition, LHCb recently measured Ar = (—0.59 ± 0.59 ± 0.21) % [s]. However, to exploit all 
available information, we use as input for the indirect CP asymmetry the result of a global 



Table 1: Experimental data on individual CP asymmetries in units of 10 



Channel 



Acp(%) 



References 



D'> -^ K+K- 



n ' n 



D° -^ K+K' 



n ' n 



D^ -^ KgKs 

D" -^ TrOyrO 

D+ -^ K+Ks 



0.00 ±0.34 ±0.13 

-0.24 ±0.52 ±0.22 

-0.43 ±0.30 ±0.11 

0.43 ±0.52 ±0.12 

-23 ± 19 

0±5 
-0.1 ±0.6 
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11 
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Table 2: Experimental averages on BR's from ref. 16 



Channel 


BR 


References 


D+ -^ tt+ttO 


(1.19 ±0.06) X 10-3 


17 


19 


D^ -^ TT + TT- 


(1.400 ±0.026) X 10-=^ 


19 


23 


^0 ^ ttOttO 


(0.80 ±0.05) X 10-3 


19 


24 


D+ -^ K+Ks 


(2.83 ±0.16) X 10-3 


E 19 


25 26 




D^ -^ K+K- 


(3.96 ±0.08) X 10-3 


19 23 


27 31 




D^ -^ KsKs 


(0.173 ±0.029) X 10-3 


J6 3 


1 33 



fit to D mixing by the UTfit Collaboration, Ay = (0.12 ± 0.12) %. The difference of the two 
asymmetries is given by 



AAcP = a'^l{K+K-) - a^^Tr+vr-) ± ^ aj?^ 



Too 



(4) 



where A{t)/Tno = {{t)K - (t)^)/rz)o = (9.83 ± 0.22 ± 0.19) % at LHCb fl]. 

Very recently, the CDF collaboration reported an updated measurement of AAqp |9]: 



AAcP = AcAK^K-)-Acp{n 



+7r- 



(-0.62 ±0.21 ±0.10) 



(5) 



with A{t)/T£,o = 0.26 ± 0.01. The individual asymmetries are listed in Table [II while the 
relevant CP-averaged Branching Ratios (BR's) are reported in Table [2] 

Combining the measurements in Table [I] with the LHCb value in eq. ([I]), with the CDF 
one in eq. ([s]) and with Ap we obtain the following average for the CP asymmetries: 



dir , 



acp(^ TT-) = (0.45 ± 0.26) % , a^l,{K+K-) = (-0.21 ± 0.24) 



(6) 
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-0.66 ±0.16) %. 



U-spin would predict acp(vr+7r^) = —aQp{K^K^). We will comment on SU{3) breaking in 
the following. 

For direct CP violation to occur, two terms with different weak and strong phases should 
contribute to the decay amplitude. For singly Cabibbo suppressed D decays such as Z^ — > 
vrvr and D — )■ KK, the CP-violating part of the relevant weak Hamiltonian is numerically 
suppressed by the ratio rcKM = ^^Vcb^ub) / (y*^Vud) ~ 6.4x10"^. Due to this suppression, the 
contribution of penguin operators is totally negligible. The possibility of direct CP violation 
then mainly rests on penguin contractions of current-current operators, which may be large 
due to Final State Interactions (FSI). Unfortunately, these long-distance effects are essentially 
uncalculable, making a prediction of a^p in these channels a formidable task. Previous efforts 



in this direction, both before and after the LHCb results, used either SU{3) [34 51 or (QCD) 
factorization 51 58 to predict a^p, or simply studied CP asymmetries as a function of the 



size of penguin matrix elements [49,59 . We improve on previous analyses in several aspects. 
First, we do not assume SU{3) nor any kind of factorization, since SU{3) appears to be badly 
broken in the decays at hand and since factorization holds only in the rric — ?■ oo limit, while 
for realistic values of rric power corrections cannot be neglected (nor estimated). Second, we 
implement unitarity constraints in a consistent way, using the wealth of experimental data on 



7iN scattering accumulated in the seventies 60 62 , yielding information on vrvr — )■ 7r7r,KK 



rescattering at energies close to the D mass scale. Third, we exploit the information coming 
from BR's to estimate the size of penguin contractions and other subleading contributions 
to the decay amplitude. Combining all information, we provide a detailed study of the 
compatibility of the Standard Model (SM) with the experimental data on CP violation. 
We conclude that, with present errors, the observed asymmetries are marginally compatible 
with the SM. Should the present central value be confirmed with smaller errors, it would 
require a factor of six (or larger) enhancement of the penguin amplitude with respect to all 
other topologies, well beyond our theoretical expectations. Thus, improving the experimental 
accuracy could lead to an indirect signal of new physics. 

The paper is organized as follows. In Sec. |2]we report the expression of the relevant decay 
amplitudes in terms of isospin reduced matrix elements and in terms of renormalization group 
invariant parameters, and give a dictionary between the two parametrizations. In Sec. |3]we 
discuss the available information on rescattering and the way to implement this knowledge 
in Z) — )■ vrvr and D — )■ KK decays. In Sec. |4] we discuss the implications of the measured 
BR's on the CP-conserving part of the amplitudes, and extrapolate this information to the 
CP-violating contributions. In Sec. [5] we present our main results on the CP asymmetries, 
and discuss theoretical uncertainties. Finally, in Sec. |6]we summarize our findings. 



2 Isospin decomposition and parameterization of D — > tttt and D — > 
KK decays 

In this Section, we write down the relevant decay amphtudes both in terms of isospin reduced 
matrix elements and in terms of renormalization group invariant (RGI) parameters, and 
discuss the relation between the two parameterizations. The isospin parameterization will 
prove useful to exploit the experimental information on final state interactions from vrA^ 
scattering, while the RGI parameterization will allow us to give a dynamical interpretation 
to the results. Before dwelling into the analysis, we give a brief summary of the relevant 
literature. 

have been used to calculate the 



Factorization approaches, such as the BSW model 52 



decay amplitudes of D decays. The experimental data favor ^ = l/N^^ ^ and demand 
significant FSI effects. Two-body hadronic D decays have also been analyzed in the diagram- 



matic approach with SU{3) flavor symmetry in refs. 38-47 , while earlier studies can be found 



in refs. 34-37 . The global fits to experimental data suggest that the color-suppressed tree is 
comparable to the color-allowed tree in size with a large relative strong phase, the exchange 
amplitude is sizable with a large strong phase relative to the color- allowed tree, and signifi- 
cant SU{3) breaking effects are required in the exchange amplitude. It is expected that the 
large exchange contribution originates from FSI. 

FSI effects on D decays have been considered in several ways: elastic and inelastic scat- 
terings, resonance contributions, etc, for example, in refs. 40,53 57,63 71 . In refs. 53 57 



Buccella et al. studied Cabibbo-allowed and Cabibbo-suppressed D decays based on a mod- 
ified factorization approximation, in which the effective parameter ^ and annihilation and 
exchange contributions are fixed from the data, and rescattering effects are assumed to be 
dominated by resonant contributions. From global analyses of the data, they showed the 
significance of the annihilation and exchange contributions and large S'[/(3) violation, where 
the latter could be explained by the rescattering effects 56 . In ref. [70], Lai and Yang con- 
sidered elastic SU{?>) rescattering (see also [72]) together with the QCDF approach for the 
short-distance annihilation amplitudes. Moreover, in refs. 63 67 , coupled- channel analyses 



of the TTTT and KK scatterings were considered for the FSI's in the D -^ tttt and KK decays. 
In ref. 



73 



Golden and Grinstein pointed out that an enhancement of CP violation in D 
decays may occur due to an enhancement of the penguin-contraction contribution as in the 
case of the AJ = 1/2 rule in kaon decays, where the AJ = 1/2 contribution dominates over 
the A/ = 3/2 one. One should note, however, that the D ^ mi data show no enhancement 
of the A/ = 1/2 over the A/ = 3/2 amplitude. We will return to this point in detail below. 



2.1 Isospin decomposition 

The effective Hamiltonian for the Cabibbo-suppressed decays with AC 
can be decomposed into A/ = 1/2 and 3/2 components, where the A/ 
involves both the current-current and penguin operators, while the A/ 



1 and AS = 
1/2 component 
3/2 component 



involves only the current-current operator 0+ = [{dL'y^lCL){uL'y^dL) + {uL'yf_iCL){dL'J^dL)]/'2. 
Namely, the A/ = 3/2 contribution involves the CKM factor V*ijVud- Denoting the isospin 
reduced matrix elements of the CP-even (CP-odd) part of the weak Hamiltonian by A{B), 
and using the original KM phase choice in which V*^Vud is real, we write the decay amplitudes 
as follows: 



A{D^ ^ n\') 






v^ 



A(D° , ;^o^on ^ ->4i3 + -^n + A + ^•'^CKM^ll + ^•'^CKM^o 
The CP-conjugate amplitudes are obtained flipping the sign of the B terms in eq. ([T]). 
2.2 Renormalization-group invariant parameterization 



In ref. 74 , a general and complete parameterization of two-body non-leptonic i?-decay am- 
plitudes was introduced based on the OPE in the weak effective Hamiltonian and on Wick 
contractions. The parameterization is independent of renormalization scale and scheme, and 
allows us to make phenomenological analyses including long-distance contributions unam- 
biguously. We apply it to the amplitudes of the D ^ mi and D — )• KK decays: 

A(D^ ^ vt+tt") = -^ [i^i(vr) + ^^(vr)] , (8) 

A{D' ^ TT+TT-) = -A, [Ei(7r) + A^in) - Pf ^^(tt) - P^'^'irc)] + A, [Pi(7r) + PgW] , 
A(D' ^ n\') = -A, [E2in) - A^in) + Pf ^^(Tr) + P^'^^in)] - A, [Pi(7r) + P^in)] , 
A{D+ ^ K+K') = A, [E,iK) - A^{K) + P,^'^{K)] + A, [E,{K) + P,{K)] , 
A{D' ^ K+K-) = A, [E,iK) + A^is, q, s, K) + P^'^{K) + P^'^{K)\ 

+ A, [E^{K) + A2(s, g, s, K) + Pi(K) + P3(K)] , 
A{D' -> ir°ir°) = -A, [A,{s, q, s, K) - A2(g, s, g, iT) + P^'^^iK) 
-X,[A2{s,q,s,K) + P^{K)] , 



where A^ = V*qVuq for q = d,b. From ref. 74,75 we have the following counting in l/Nc. 



El and Ai are the leading amplitudes, all other amphtudes are suppressed by 1/Nc except 



for P3 and P^^^, which are suppressed by l/N^- The amphtude for D^ — > K^K^ is 1/Nc 
suppressed and originates from S'f/(3) breaking effects. 

In terms of the RGI amphtudes, neglecting the contribution proportional to tckm to the 
A terms, the isospin amplitudes can be written as 
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Ad[^i(7r) + £;2(vr)], 



A, [2Ei(7r) - ^2(vr) + 3A2(7r) - W^^^^^-k) - 3P^'^'{n 



(9) 
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V3 

-V3X4Piin) + Psin)], 
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\K 



•^11 — ^a 



E,{K) - '^-A,{K) + \A,{q, s, g, K) + P^'^{K) 



Bf, = \4E,iK) + P,iK)] 



A^ 





B, 



-A, [E,iK) - A,iq, s, q, K) + 2^2(5, g, s, K) + P^'^^^K) + 2P^''^{K)\ 
-X4E,{K) + 2A2{s,q,s,K) + P,{K) + 2P3{K)]. 



Therefore, one expects Bq to be 1/A^c- suppressed with respect to Aq. This suppression is 
partially compensated by the Clebsch-Gordan coefficients, so that the two amplitudes could 
be of the same size. Concerning the amplitudes with kaons in the final state, they are all 
leading in the l/N^ counting; however, a cancellation between the emission and annihilation 
parameters may occur in Ai^, possibly leading to an effective l/N^ suppression. 

Neglecting the 0{1/N^) contributions, the combinations of the effective amplitudes for 
the TTTT modes are written in terms of the isospin amplitudes as 



(10) 



E,{7r) + E,{n) = -X-\ -Ah 



,{n) + A,{n) - Pr\^) = 


I / A^ \ 


2{7r) ~ A2{n) + P^'^\n) = 


-\4^[V2a; + a-, 


Pii^) = 





while those for the KK modes are given by 



A,{K) = \4 \ (-2 <3 - A^, + Sfi - < + B^) 
A2(g, ., g, K) = \4 \ (-<3 + <, -B^,+ < - B^) , 
A,{s,q,s,K) = \4\{-Bf,-B^), 



;ii) 



6 



E,{K) + P,^'^iK) = Xf \ (-^f; + A^, + B^,-A^ + B^) , 

Before turning to the phenomenological analysis, we discuss the constraints imphed by 
unitarity on the isospin aniphtudes. 

3 Rescattering and unitarity 

Unitarity of the S-matrix imphes constraints on weak decay matrix elements, provided that 
the strong 5* matrix at the relevant energy scale is experimentally accessible. As we discuss 
below, this is indeed the case for Z) — t- tttt and KK decays, leading to interesting constraints 
on the decay amplitudes. 



Notice that any Wick contraction, as defined in refs. 74 , 76 , can be seen as an emission 



followed by rescattering 76,77 . Thus, rescattering establishes a link between emissions and 



long-distance contributions to other subleading topologies such as penguins. 

3.1 Coupled-channel unitarity 

We split the effective Hamiltonian for weak charm decays into a CP-even H^ and a CP-odd 
Hi part. Then we can write 



The S matrix can be written as 



(12) 
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( D^D 

TTTT — )• D 

KK ^ D 

V 



D ^7171 D ^ KK 

TTTT — )■ TTTT TTTT — )• KK 

KK^nn KK -^ KK 



...\ 



/ 



-zCP(T) 



Ss 



(13) 



where the time reversal of T is equal to the CP conjugate of T: T(T) = CP(T) = T^ — iT^, 
and 5*5 is the strong interaction rescattering matrix. Unitarity of S (and of Ss) implies, at 
lowest order in weak interactions, 



T 



R 



Ss{T 



R\* 



T^ = SsiT 



i\* 



(14) 



where separate equalities hold for T^ and T^. These equalities can be used to reduce the 
number of unknown hadronic parameters in the decay amplitudes, if Ss is known indepen- 
dently. The simplest case corresponds to decay channels where Ss can be approximated with 
a pure phase e^*''. Then we obtain 



T 



R 



\T^\e'\ 



T' 



i\js 



T' e 



(15) 



where 5 + vr is also possible. Notice that, even in this simple case, eq. (15) cannot be used 
to add FSI to factorized amplitudes, since the identification of factorized results with \T\ (or 
Re T) is ambiguous. 

In principle, this could be the case for 1 = 2 S-wave tttt — t- tttt scattering, whose phase 
can be extracted from the data in ref. ||62|: 

C' = (-8 ± 5)° . (16) 

However, there is a sizable inelasticity in this channel at the D mass, so that the information 
above cannot be used (see the discussion below on multi-channel unitarity). This does not 
spoil the rescattering analysis since, as we show below, the D ^ tttt BR's fix the relative 
phase of A2 and ^q with an excellent accuracy. 

Concerning the 1 = 1 KK rescattering, if it were elastic we would have arg^{| = 
Big All = argi3{^ up to a vr ambiguity, leading to the absence of direct CP violation in 
D^ — )■ K^Ks- However, it is well conceivable that KK scattering at the D mass is inelastic, 
so that we do not impose the relation above. 

The case of / = amplitudes is more involved. Experimental data on vrvr and KK final 



states have been collected in refs. 60 and 61 respectively. The data on vrp — )■ KsKgn show 



a strong suppression of the vrvr — >■ KK amplitude at energies close to the D mass, as can be 



seen for example in Fig. 6 of ref. 61 . Conversely, the extraction of isospin amplitudes from 



data on tttt — > tttt scattering at the D mass is ambiguous, leading to widely different results 



for the inelasticity. For example, ref. 78 provides four different amplitude fits corresponding 
to discrete ambiguities; one of them gives results compatible with the KK data close to the 
D mass, while the others point to violations of two-channel unitarity. The latter could be 



due to the four pion channel, see for example Fig. 3 of ref. 79 . Thus, two scenarios may be 
envisaged. 

3.1.1 Two-channel analysis of J = amplitudes 

First, one can assume that the strong 5* matrix is well described by a two-channel analysis 
with TTTT and KK states only. Indeed, two-channel fits give a reasonable description of data in 
a wide range of energies (see for example Fig. 1 of ref. (SOl). The corresponding two- by-two 
symmetric rescattering matrix can be parameterized as 



where r) is the inelasticity parameter. We extract the 1 = S-wave scattering phases of tttt — )■ 
TTTT and TTTT — )■ KK and the inelasticity parameter rj at the D mass from the experimental 



data in refs. 60 61]: 

5i = (40±10)°, (51 + ^2 = (360 ±60)°, r/ = 0.95 ± 0.05 



where the last value has been estimated using KK data. 

Unitarity implies the following equation for / = CP-even amplitudes: 



A^ 



Tje 
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^2 gJ(i5i+<52) 
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2i&2 




■K,K 



and an identical equation holds for the CP-odd amplitudes i3g' . Defining ^q 



and A^ = 1^0 1^*"^° ) 61- (19) can be written as 



sin(v9^ - ^2) 
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1 — rj 
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sin(<^o - K 



cos(v9o - Si] 



(19) 

\A^\e'^o 

(20) 



where one can add vr to the phases ip^ and (p^ simultaneously. From these equations, we find 
that the ratio |^^/^q| obeys the following constraints: 



1 — ?7 
1+V 



< 



K 



A^ 
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< 



l + ri 



1 



T] 



(21) 



and the phase differences (p^ — 5i and ip^ — 82 are determined in terms of |v4o /v4o| and r] 



In the limit of r] — )■ 1, where the scatterings are elastic, eq. (|19|) can be written as 



e-"i \Al\e 
and the strong phases are then given by 

ifl = 5i + mr 



\A, 



V^ 



KlJip^ 



^2i&2 \ aK\ ^-i<fio 



\A^ 



62 + mvr . 



(22) 



(23) 



where n and m are arbitrary integers. Similarly, the strong phases of the CP-odd amplitudes 
are given by (pQ = 61 + n'lT and (pQ = 62 + mV, where n' and m' could be different from n 
and m. In this case, CP violation cannot be generated from the interference of Aq [Aq) and 
Bq (Bq). Thus, in this scenario, given the small inelasticity of tttt scattering, we expect that 
CP violation in D — )■ vrvr decays mainly arises through the interference of Bq with A^- 

3.1.2 Three-channel unitarity 

In the second scenario, instead, we allow for a third (effective) channel to give a sizable 
contribution, thus reconciling the large inelasticity solutions of vrvr — > vrvr amplitude fits with 
the KK data. This corresponds to a three by three Ss matrix in which the KK channel 
is almost decoupled, leading to a situation similar to the one described above but with vrvr 
coupled to the third effective channel with a large inelasticity (small r]). If the KK channel 
is decoupled, unitarity fixes the phase of Aq and Bq to be equal to 62 + m^ j8l]. Conversely, 
since the vrvr channel has a large inelasticity, the solutions of two-channel unitarity discussed 
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lA^I [GeV] 



|A;;| [GeV] 



arg(A^/A5 [-] 



Figure 1: From left to right, p.d.f. for 1^21' l-^ol ^^^ ^''^§,{•^2 / -^o) i^ ^^^ three-channel sce- 
nario. Here and in the following, darker (lighter) areas correspond to 68% (95%) probability 
ranges. 



above give essentially no constraint on absolute value and phase of Aq and Bq. Thus, in this 
case CP violation in Z) — )■ tttt can also arise from interference between Aq and Bq. 

We have checked numerically that the results obtained in the three-channel scenario are 
essentially identical to the ones obtained in the most general case, where more than three 
channels contribute to the rescattering so that no significant constraint can be obtained from 
unitarity. 

4 Branching ratios and CP-even contributions 



The decay width of the process D — t- PP is given by 

Pc 



T{D^ PP) 



87rm|) 



\A{D^ PP)\' 



(24) 



where pc = A/m|^ — 4m|^/2 is the center-of-mass momentum of the mesons in the final state, 
and an extra factor 1/2 must be added in the case of D^ — )■ 7r°7r°. We adopt rriK = 0.498 
GeV, m^ = 0.135 GeV, m^ = 1-865 GeV, r^o = 410.1 x lO^^^ sec and r^± = 1040 x 10"^^ sec 
in numerical analyses. In this Section, we discuss the determination of CP-even amplitude 
parameters from the measured BR's reported in Table [2j Here and in the following, we 
follow the inferential framework outlined in ref. 82 . In particular, we obtain the 68% and 



95% probability regions by integrating the posterior p.d.f. around the most probable value(s). 
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4.1 7171 isospin amplitudes 

In the case of D — )■ vrvr, the BR's are sufficient to determine |^o 2I ^^^ ^^e relative phase. The 
magnitude of the 1 = 2 CP-even tttt amplitude A2 can be extracted from BR{D^ — )■ vr^vr^): 



and then ^q and the relative phase can be obtained from BR{D^ -^ vr^vr'^) and BR{D^ — )■ 
7r+7r~). From the probability density function (p.d.f.) in Fig. Q^we obtain 

l^^l = (3.08 ± 0.08) X 10"^ GeV , (26) 

l^^l = (7.6 ±0.1) X 10"^ GeV, 
arg(^^M^) = (±93±3)°. 



Notice that the results in eq. (26) exclude order-of- magnitude enhancements of the / = 



amplitude. The quality of the fit to the BR's is excellent. 

4.2 KK isospin amplitudes 

In the case of D — ?► KK decays, the BR's are not sufficient to determine all isospin amplitudes. 
Given a value of Aq that satisfies the unitarity constraints, we solve for |^^/2 + A^^l, 
\Af^ - A^sl and arg((^fi - ^{|)/^^) using the three BR's. The p.d.f. for \A^J2 + ^fj, 
|^{^— ^^1 vs l^^l and aig{{Aii — Af^) / A^ ) are reported in Fig. |2} In order to reproduce the 
CP asymmetries, the degeneracy in arg((^{^ — A^^)/Aq) is broken, with a mild preference 
for the negative solution. 

An interesting result is given by the CP-conserving contribution to BR{D^ — )■ K^K^), 
which should vanish in the SU{?)) limit. We obtain instead a result comparable to all other 
amplitudes in the KK channels (see Fig. [2]): 

I^II -•^ii--K\ = (5-0 ± 0.4) X 10^^ GeV , (27) 

showing explicitly a breaking of 0{1) of the S'f/(3) flavour symmetry. Also in this case, we 
obtain an excellent ffi of the BR's. 

4.3 RGI parameters for CP conserving contributions 



From eqs. (10) we obtain the following results for the pion RGI parameters in the three- 



channel scenario: 



Ei{7i) + E2{7x) = (1.72 ± 0.04) X lO^^ e'^ GeV , (2^ 

Ei{7t) + A2(7r) - Pf i^(7r) = (2.10 ± 0.02) x 10"^ e^{5±(7i±3)°) Q^y ^ 



11 




0.2 0.4 0.6 0.8 1 1.2 1.4 



xlO 



lA^'l [GeV] 








xlfl' 










>, 














.■^ 














w 
















20 


- 










o 














>. 














.■ti 














O 


15 


- 










(Tl 










I 




J3 










I 




O 










■ 


















CL 


10 

5 


- 


1 


1 


1 






) 


0.2 


0.4 


0.6 


0.8 1 



>> 10000 



|A^^/2+A^,| [GeV] 




A^ +A^^+A„K| [GeV] 



Figure 2: From left to right and from top to bottom, p.d.f. for |^{\— ^{|| vs \Aq\, arg((^{\ — 
Afr^)/A^), |^{|/2 + Afi\ and \ - Af^ + A^i+ A^\ in tlie two-cliannel scenario. In tlie tliree- 
channel scenario one obtains essentially identical results. 



^2(vr) - A2' 



vr 



+ Pf^^(7r) = 



= (2.25 ±0.07) X io-6e*(5T(62±2)°) Q^y 



with a two-fold ambiguity and generic 6. These results show that the Ei (it) parameter does 
not dominate the decay amplitude, and that 1/A^c-suppressed topologies are comparable to 
Ei{7[) with a large strong phase difference (this is evident by comparing the second and third 
lines of eq. ([28|)). This also shows that power-suppressed amplitudes in the mc —)■ 00 limit 
are of the same size of leading ones. 

Let us now turn to the KK channels. The result in eq. (27) implies, using eq. (tsl), that 
the S'f/(3)-suppressed combination of subleading amplitudes ^2(3, q, s, K) — A2(g, s, g, K) + 
pGiM^j^^ is of the same order of the leading contribution Ei[K). 

We conclude from the analysis of D — )■ vrvr and D — > KK BR's that subleading topologies 
are of the same order of leading ones, with a breaking of S'f/(3) of 0{1). This is the starting 
point for our study of CP-violating asymmetries in the next Section. 
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Figure 3: From left to right, p.d.f. for the prediction (first row) and fit (second row) of Aa^p, 
a^p(7r+7r^) and aQp{K^ K~) in the two-channel scenario. Third and fourth rows: same as 
the first two rows in the three-channel scenario. Darker (lighter) areas correspond to 68% 
(95%) probability ranges. The dotted (dashed) lines correspond to 68% (95%) experimental 
ranges from eq. ^. 



5 CP asymmetries 

We turn to the main point of this work, namely the attempt to estimate the possible size 
of CP asymmetries in the SM and to quantify the agreement of the SM with experimental 
data. 

Before dwelling in the analysis, we remark a few relevant points: 

• present experimental data point to a larger CP asymmetry in the tt+tt^ channel with 
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Figure 4: From left to right, p.d.f. for the fitted Aa^p for different values of k in the two- 
and three-channel scenario. All the p.d.f.'s have been scaled to fit in the same plot. 



respect to the K^K one (indeed, the latter is compatible with zero at less than la); 

CP violation is always proportional to subleading contributions; in the case at hand, 
CP asymmetries in the K^K~ and vr+vr" channels are due to penguin contractions of 
current-current operators, while in K^K^ also annihilations contribute; 

in the two-channel scenario, one has to a good accuracy argSJ = arg^Q and CP 
violation can occur only through the interference of Bq with A21 leading to a suppression 
of the CP asymmetry with respect to the three-channel scenario; 



given our phase convention for the CKM matrix, CP violation in the tt+tt channel is 
signaled hy Bq ^ Q, while in the K^K~ channel one must have B^^ ^ Afi — Ai^ or 

K 

■ 



B^ ^ Af 



Thus, to estimate CP asymmetries we need to estimate the size of subleading amplitudes. 
From the analysis of the BR's presented above, we do not see any evident suppression of 
subleading terms, so that we impose generically 



\Bl\<K\Al\, 
\B^ - A^\ < k\A^\ 



(29) 



\B^ 



\AyY Ayj,)\ < l^\Aii Ai^\ , 



\Ki 



where k parameterizes the size of the subleading terms. In terms of RGI parameters, this 
amounts to 



\Pi{n)\<K 



Ie^{7c) - ^E^in) 



A,{n) - Pl^'^in 



(30) 



GIM/ 



Pi{K) - P^^^(K) + A^iq, s, q,K)\< k\E,{K) - ^^(g, s, g, K) + 2^2(5, g, s, K) + P^^^(K)| 
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Figure 5: First row: p.d.f. for the parameters in eq. (28). Second row: p.d.f. for |Pi(7r)| in 
the two- and three-channel scenario. All the p.d.f. 's have been scaled to fit in the same plot. 



\P^iK) - Pr\K) - A^iq, s, q, K)\ < k\E,{K) + A,{q, s, q, K) + P\^'^{K)\ , 

where P3 and P^^^ have been neglected. 

Let us first present the results for 1 < k < 8 and then comment on the values of k that 
we consider acceptable. We can follow two different avenues. The first possibility is to give a 
prediction of the CP asymmetries as a function of k and compare it with experimental data. 
The second option is to fit the measured CP asymmetries as a function of k. In this case we 
can also study the values of the subleading topologies selected by the fit and compare them 
with our (albeit vague) theoretical expectations. 



In the upper part of Fig. pi we present the predictions and fit results for Aa^p 



IqP 



TT ' TT 



and a'^p{K^ K~) in the two-channel scenario. We see that the generic prediction would give 
much smaller asymmetries, and that the prediction does not reach the present experimental 
value within 2cr for values of k < 8. In the three-channel scenario, instead, we obtain the 
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Figure 6: P.d.f. for the CP asymnietries in the two-channel scenario for different values of k. 
All the p.d.f.'s have been scaled to fit in the same plot. 

results in the lower part of Fig. [3j Since in this case the pion amplitudes are less constrained 
by unitarity, the predicted asymmetries are larger than in the two-channel scenario, and the 
present experimental value can be reached within 2a for k > 5, but even for k = 8 the 
prediction is still la from the experimental result. The p.d.f. for Aa^p for different values 
of K, can be found in Fig. |4j 

To assess the compatibility of the experimental result with the SM, we can compare the 
distribution for -Pi(vr) obtained from the fit for different values of k with the distribution of 



the pion amplitude parameters obtained from the BR's in eq. (28). To this aim, we report 
in Fig. Is] the p.d.f. for the absolute values of the parameters in eq. (28) and for Pi(7r) 



for different values of n in the two scenarios. We notice that in the three-channel scenario 
the preferred value for |Pi(7r)|, corresponding to the central value of the measured Aa^p, is 
around 1.3 x 10^^ GeV, about 6 times larger than the RGI parameter combinations obtained 
from the BR's. In the two-channel scenario, instead, even for k = 8 the fit is still pulling 
|Pi(7r)| to the upper edge of the allowed range, showing that the present central value cannot 
be reasonably accommodated in this scenario. 

For the sake of completeness, we report in Figs. [6] and [7] the p.d.f.'s for the fitted CP 
asymmetries for different values of k. 
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Figure 7: P.d.f. for the CP asymmetries in the three-channel scenario for different values of 
K. All the p.d.f.'s have been scaled to fit in the same plot. 

6 Summary 

We have analyzed the D — )■ KK and D — t- tttt decays within the SM, assuming only isospin 
and using the information from tttt scattering and unitarity. We have considered two possible 
scenarios for the strong S matrix (two- and three-channel unitarity). We have performed a 
fit of the CP conserving contributions from the CP-averaged BR's, obtaining information on 
isospin amplitudes and RGl parameters. We have predicted and fitted the CP asymmetries 
in the two scenarios. 

Considering the more conservative three-channel scenario, we conclude that, with present 
errors, the observed asymmetries are marginally compatible with the SM. This conclusion 
holds also for the most general scenario with even more coupled channels in the 1 = 
rescattering, where no significant constraints arise from unitarity. Should the present central 
value be confirmed with smaller errors, it would require a factor of six (or larger) enhancement 
of the penguin amplitude with respect to all other topologies, well beyond our theoretical 
expectations. Thus, improving the experimental accuracy could lead to an indirect signal of 
new physics. 
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